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Abstract 

We consider perturbations in the fields of the 't Hooft-Polyakov monopole, dyon 
and point electric solutions. We find a series of bound modes where the fields are 
confined to the core of the respective solutions, and these are interpreted as bound 
states of a gauge boson with the respective solutions. We discuss the spectra of 
these bound states of the various two particle systems of the Yang-Mills-Higgs 
model in connection with the duality of Sen. 
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1 Introduction 



In Grand Unified Theories with a simple compact gauge group the series of symmetry 
breakings necessary to reach the standard model inevitably gives rise to monopoles. 
These monopoles interact with the other particles of the theory and the nature of this 
interaction is of interest. This interaction has been studied by a number of authors 
10]. As well as scattering states bound states of bosons and fermions with monopoles 
have also been found Perhaps the most interesting finding is the baryon decay 

catalysis by a monopole 

Here we are most interested in the work of Sonoda In ref. [§] the physics of 
the massive gauge boson (W^) and the neutral Higgs scalar (H) in the presence of a 't 
Hooft-Polyakov monopole p, 3 was investigated. Spherically symmetric perturbations 
of the fields of the 't Hooft-Polyakov monopole were considered and a series of bound 
states were found. The eigenvalue of the lowest bound state was numerically estimated 
to be ^ 0.876m in the Prasad-Sommerfield limit (m is the mass of the gauge boson). 
We investigate these bound states further by using a more systematic approach to the 
study of general angular momentum perturbations in the fields of the 't Hooft-Polyakov 
monopole. Our motivation for studying them is, however, quite different. 



In ref. [10| Montonen and Olive observed that in the Bogomol'nyi limit [11| the 



masses of all the particles of the theory are given by the same universal mass formula 



m 



^ = rj'^iq'^ + g"^), where q is the electric and g the magnetic charge of the particle. 
This formula applies to both the classical non-singular monopole and dyon solutions as 
well as the normal particles, and in part motivated Montonen and Olive to conjecture 
that the VT-bosons and the monopoles are dual particles. The recent construction of 
supersymmetric models ||12|, [1^ which realize this electromagnetic duality has stimulated 
a great deal of interest. 

We do not consider these supersymmetric models here, however, but consider the 
non- supersymmetric model of 't Hooft and Polyakov. We note from the literature cited 
above that the massive particles of this theory, whose masses are given by the universal 
mass formula, namely the monopole, dyon and ly-boson form bound states with one 
another. So here we seek to answer the question: are the bound state energies of the 
various two particle systems related by electromagnetic duality? 

To answer this question we need to find the values of the bound state energies in 
question. The bound state energies for the monopole-dyon system were found in ref. O 



Here we try to complete the picture by using perturbative methods to investigate the 
bound state energies of l^-bosons in the presence of a monopole, a dyon and another 
VT-boson. 
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2 Yang-Mills-Higgs model 

We are considering field configurations in the Yang-Mills-Higgs theory with the lagrangian 



2\2 



where 



and and $° are SU{2) triplets. We will consider monopole, dyon and point electric 
solutions to the field equations 

It will be convenient later on if we use dimensionless fields and co-ordinates, so we make 
the replacements 

r = qr]R, 0" = 77$^ = r]A''^. 

Note that the mass of the massive gauge bosons is = qt]. With these dimensionless 
variables we will take the ansatz for the monopole with unit magnetic charge to be |[ 



r 



The profile functions are solutions to 



r'^ dr ^ dr^ ~^ ~'~ 2 ^ 

d^W o W(W-2) 

' H\W -l) + iW -l) ^ ' = 0, 



dr'^ 

where we have defined j3 = 2\/q^. We will only be considering the Bogomol'nyi limit 



(/3 — > 0) where we have the analytic Prasad-Sommerfield solution |T5 



sinh(r) r' sinh(r) 
For the dyon with unit magnetic charge we will take the ansatz 



r 
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where the profile functions are solutions to 



ar \ ar r^ I 



+ iH^-J^)iW-l) + iW-l)^^^^ 







ar \ ar r^ 



In the Bogomol'nyi hmit we will follow Bogomol'nyi |TT| and take the dyon solution to 
have the form 

0« = -C2 x), W^ = ^/l- C2 W"(Vl - C2 x), PVo'' = C0"(Vl-C2 x) 

where (0^, W^^) is the Prasad-Sommerfield monopole solution given above. The parameter 
C is given by 

C = , where e = - — 

and Q is the electric charge of the dyon. The mass of the dyon in dimensionfuU units is 
given by 

AT^rj 



Q 

and setting e = gives the mass of the monopole. From this we can see that the quantity 
(1 — C2)^/2 = (^]^ _ g2~)-i/2 appearing in the rescaling above is the ratio of the monopole 
and dyon masses. Note that with this form of the solution the profile function J is given 
by J = CH. 

To investigate the bound state energy of a IV-boson in the presence of another W- 
boson we will consider perturbations about the point electric solution 

^' = {i-h n = --^- 

Substituting this ansatz into the field equations gives 

D\(f = AnZ6^{r) 
DiGl = ATiZ5\r) 

and so we have a point electric source of strength AttZ in rescaled units. Note that in the 
rescaled units the VT-boson has charge 1, and so Z = l/47r corresponds to a W^-boson 
source. 

Now to find the equation for the perturbations (50", SW^) about a solution (0", W^) 
to the field equations, we take the second functional derivative of the action. To remove 
the gauge degrees of freedom we use the background gauge condition 

D^SW^ + e"^'=0''50'= = 0. 



Imposing this gauge condition, we find the gauge fixed perturbation equations are 

{D^)ab + 5a6(f (0 ■ - 77') + ■ 0) + (/? - 1)0'^0' 2e'^^'=Z}^0^ \ ( ^"^^ \ (U 

= 0, 

where the metric is given by g^'^ = diag{— 1,1,1,1). The gauge condition does not, 
however, fix the gauge completely. To remove the remaining degrees of freedom requires 
the introduction of Fadeev- Popov ghosts rj"" which satisfy 

{{D'),, + 6ab<P ■ - 0'^0') v' = 0. (2) 

If we count the number of degrees of freedom we find that equation |1] has 15 and equation 
0has 3. When calculating any quantum corrections the Fadeev-Popov ghost contribution 
appears with a factor of —2 relative to the contribution from the modes of equation |I[ 
In this way the Fadeev-Popov ghosts cancel 6 degrees of freedom of equation |l| leaving 9 
physical degrees of freedom. They correspond to 6 degrees of freedom for the 2 massive 
gauge bosons, 2 degrees of freedom for the massless gauge boson and one degree of 
freedom for the Higgs boson. 

We seek to solve equation ^ for the physical bound states for general angular mo- 
mentum in the presence of the three solutions given above. To check that the states we 
find are actually physical states we shall also solve equation ^ If the eigenspectrums of 
equation | correspond to an eigenspectrum of equation |I] then that eigenspectrum is not 
physical; it is cancelled by the ghosts. 



3 Perturbations about the monopole 

If we substitute the monopole solution into the gauge fixed second variation of the action, 
we find 



ss, 



9f 



- I r'^dr / dn 



Dlb + 5"' (^iH^-l) + hA + f''f\/3 - 1)H 



2/ xab -a-6^ 



2/xab _ f.o.f.b\ 



(3) 



and the Fadeev-Popov ghost term is 



5S 



FP 



r'dr / d^T]" Di^ + H'{6' 



2/xab _ fa^b\ 



(4) 



where 



Sabdo + Sab 



w 



-V' + ^ + 1 - 2—e,abLi + f'^f^— 
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— h -^{w - 2) - — 



(DicPY = r'r^(H' -—(l-W)] +6'''—(l-W) 
\ r J r 



and L is the angular momentum operator L = —if x V. Note that because the monopole 
solution is time independent and Wq = we can separate out the time derivatives and 
set up a standard eigenvalue problem. 

The monopole solution is invariant under rotations generated by = L + T + S*, and 
so we seek to expand the perturbations as eigenstates of K. First we consider the scalar 
fields {Wq, 77", 0") which have spin zero and so K = L + T. Simultaneous eigenfunctions 
of K^, K3, L^, and can be constructed from eigenfunctions of L^, L3, and T3 by 
using standard Clebsch-Gordan methodology 

= E C!ZimyimYtm, where m = rrii + rrit 

mimt 

where C are Clebsch-Gordan coefficients and 1^^, are the spherical harmonics. We will, 
however follow ref. []T7| and define vector spherical harmonics in terms of operators acting 
on the spherical harmonics by 

^ kjm k,j '~^j^k,m- 

The vector operators Oj are given by 

00 = iL 

01 = f{L + l)+ifxL 
O-i = fL — if X L 

where 

/ 1 



4. 

The normalization constants Nkj are given by the orthogonality condition 



The vector spherical harmonics Y^j^ defined above are eigenfunctions of K^, K^, L^, and 

^^3^ kjm ~ '"'^kjm 

L'Y,)m = (k + j)ik + J + 1)Y,)^ 

kjm kjm. 
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We now expand the scalar fields as 

km i=— 1 

E E 

km J=— 1 

E E 

km j=— 1 

(5) 

Note that the orbital angular momentum I is given hy I — k -\- and so for each value 
of k there will be three values of So the subscript j on the vector spherical harmonics 
denotes the three different values of / that are possible for each value of k. 

For the spatial components of the gauge field we have non-zero spin and isospin and 
so K — L-\-T-\-S. There are two different ways of adding three angular momenta, 
differing in the order in which the addition is performed. We choose to first construct 
eigenvectors of 

J^T + S. 

Note that J here is not the total angular momentum operator. For spin one and isospin 
one there are three possible values of j, j = 0, 1,2 denoting scalar, pseudovector and 
tensor representations of J. If we then consider eigenstates of 

K ^ J + L^T + S + L 

we can see that in general for each k\ for j = there is only one value of /, / = A;; for 
j — 1 there are three values of/, / =|A; — 1|,A;,A;-|-1; and for j — 2 there are five values 
of Z, I — |/c — 2|, |A; — 1|, /c. A; + 1, /c + 2. So we introduce the variables ji and j2 where ji is 
the J eigenvalue and j2 takes the values j2 = — ji, ■ ■ ■ ,ji so that in general, I is given by 

1 = k + j2. 

With this notation we expand in tensor spherical harmonics as 

2 ii 

km 31=0 32=— 31 

where the tensor spherical harmonics Y^^^^^^ are defined in terms of operators acting on 
spherical harmonics by 

yai _ TSJ-'^l'^rjai y 
k3ij2m ^ " kjij2 ^3132 ^ ' 
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where the tensor operators Oj^^j^ are constructed from the Oj. The tensor spherical 
harmonics defined in this way are eigenf unctions of K'^, K3, and 



kj\j2m 

K Y""^ 

3 kjij2in 

J 2 yai 

kjij2rn 

j2-yai 

kjij2m 



k{k + 



kjij2m 



kjij2m 



]2m 



kjij2in- 



Note that the scalar (ji = 0), pscudovector (ji = 1) and tensor (ji = 2) representations 
of J have all been incorporated into the single set of functions l^jj^jm- 

Now the monopole is also invariant under the action of the parity operator P. If we 
apply the parity operator to the perturbations, then from 



pya 

kjm 



'^-i^^+jy^ onH PV"' — (_Y\k+j2yai 

K ^) ^ kjm '^'-^^ ^ kjij2m \ ^ kjij2m 

we find that the perturbations separate into two channels under parity given by 



km 



1 


\ 




/ 


\ 


-U/km 
















^km 








yykm 
















) 






/ 



where the superscripts are such that 

+ 



denotes parity 
denotes parity 



-1 



The perturbation equations for ip^^^ decouple. Note that the vector and tensor spherical 
harmonics are eigenstates of K^, K^, and and so are not eigenstates of the charge 
operator 

Q = r.T. 

Substituting these expansions into 5Sgf and performing the angular integrals gives the 
perturbation operator 



M, 



(±) 



a/3 



w 



(±)(1) 



+ 



(±)(2) 
:/3 



+ 2i^::cS«'' + 2^(H---2)cS'« 



r 



(±)(6) 

:/3 
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where (±) denotes the two channels with different parity, the subscripts a and /3 denote 
the modes labeled by ji and j2 (i-e. a = (ji, ^2) and (3 = {ji,j2)), and is a projection 

y^i^ are symmetric 



operator which is 1 for a (j) field and for the W fields. The Cl^^*^^ 



matrices of coefficients resulting from the integration over angular co-ordinates. This 
notation is chosen to be the same as that used in ref. [O, and the matrix elements of 



C^^^^~^^ are taken from there. 

The perturbations for the time components (and the ghosts) similarly separate under 
parity 



km 



( W^"" 



km 



where the superscripts ± denote the parities as before. 

The perturbation operator for the time components (and the ghosts) is 



M, 



(±) 



a/3 



+ 



2W 



+ 



a 



r.2 af3 



(±)(1) 



+ 



where the C, 



(±){1,2) 
a/3 



elements are the same as those for the 6 terms. 




k = modes 

If we turn our attention to the k = modes we find that, because of the constraints 
I = k + j2 > and A; > 0, the number of channels for k = are reduced. For k = the 
allowed perturbations are 

~ V ^00 / 

The ipQQ mode satisfies the perturbation equation 
where 

.2,6 ^ 5W\ 1 ^2 , 8W\ 2W{W-2) 



Mn = -Vt + --6— + -— + -H' + -— + - 2 



3 3 3 r 3 
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and the ■j/'oo mode satisfies the perturbation equation 



2 ' 

If we substitute in a numerical solution for the profiles H and W, these equations can 
be solved numerically for general (3 by standard matrix methods. The most interesting 
case is the Prasad-Sommerfield limit where we have an analytic solution. Despite the 
existence of the analytic solution we have been unable to solve the perturbation equations 
analytically. We can however, proceed analytically if we take the far field forms of the 
profile functions 

H^l--, W^l for r^oo. 
r 

Substituting this into the perturbation equation for ipQQ gives for r ^ oo 



for which we can find two orthogonal eigenvectors. One of the eigenvectors is 



where the function wi{r) satisfies 



This has spherical Bessel function solutions for all o;^ > 0. The other eigenvector is 



ws\ / 1 \ 



and the function W2{r) satisfies 



-V^ + 1 - ^2 = U}'^W2. 

If we write W2{r) — u{r)/r, rescale p — 2(1 — uP'^l'^r and define A = 1/(1 — uj'^Y^'^ then 
we obtain Whittakers form of the hypergeometric equation 



' 1 A \-p^\ , , 



with /i^ = 1/4. For c<j^ > 1 this equation has Whittaker function solutions Wa,^(p) for 
all A which are oscillatory at large p. For 1 > a;^ > the Whittaker function solutions 
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only exist for A = 1,2,3,... and the solution decays exponentially at large p. These 
bound states are analogous to the states of the non-relativistic hydrogen atom, and have 
eigenvalues 



ujI = 1 t: where n = 1, 2, 3, 



1 

The Whittaker function solutions can actually be written as Laguerre polynomials, since 
if we set u = pe'^^'^v we get Laguerre's equation 



(p^ + (2-p)^ + (A-l))t;(p) = 0, (7) 

which for A = n E Z has associated Laguerre polynomial solutions L\{p). 
If we now consider the equation for i^qq in the limit r ^ oo we find 

( -Wl + l-l \( W^oo \ 2 ( ' 








2 ( ^^11 



The Higgs scalar field has spherical Bessel function solutions at large r, while the W^f 
field has the same type of bound solutions as above with eigenvalues 

cu^ = 1 — \ where n = 1, 2, 3, . . . 

These spherical Bessel function and Whittaker function solutions are finite at the 
origin and so are physically acceptable solutions for all r. They were obtained using the 
far field limit for the profiles which ignores any effect of the monopole core, so it is not 
obvious that the eigenvalues obtained will still hold when we use the full analytic solution. 
However, numerical solution of the perturbation equations with the full analytic form of 
the profile functions gives numerical results which are consistent with the analytic result 
obtained in the far field limit. 

Table 5 shows the analytic and numerical results for the eigenvalues of the lowest 
identifiable bound modes. The numerical results were obtained by solving the V'oo equa- 
tion using standard matrix methods. The non-linear map s = tanh(r) was used to assign 
more lattice points to the core of the monopole where the fields are varying most. Apart 
from the n = 2 eigenvalue we can see that the agreement between the analytic and nu- 
merical results is fairly good. In trying to argue that the numerical results actually agree 
with the analytic result we have to address two principle shortcomings of the numerical 
results: the difference between the n = 2 numerical and analytic results; and the limited 
number of bound modes identified. 

The n = 2 bound mode differs from the other bound modes in that the fields of this 
mode vary most in exactly the same place that the monopole profiles vary most. We 
would argue that this leads to the error in this mode being greater than for the other 



modes. In fact we found a similar situation in ref. |]18[ where we considered perturbations 



about a Nielsen- Olesen vortex 19 . There we found a similar bound mode whose fields 
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varied most in exactly the same place as the string profiles, and its eigenvalue was the 
most sensitive of all the modes to the accuracy of the string profiles used. 

This n = 2 mode also appears to be the one greatest affected by the omission of 
the origin, which is necessary because of the centrifugal terms in the perturbation 
equation. The effect that the omission of the origin can have is best illustrated by 
considering the analogous case of perturbations about the kink [^. In this case the 



perturbation equation can be solved analytically and two discrete modes are found. One 
is a cj^ = translation mode and the other is a ti;^ = 0.75 bound mode (m is the mass 
of the Higgs boson in this case). If we solve this eigenproblem on a lattice numerically 
we obtain agreement with the analytic values for the discrete modes above. However, if 
we remove the lattice point at the origin, there is a 0.01 shift in the eigenvalue of 
the bound mode, whereas the = translation mode eigenvalue is unaltered to order 
0.01 m^. 

We would therefore argue that these affects, notably the omission of the origin, are 
responsible for the disagreement between the numerical and analytic results of the n = 2 
mode, and that its true eigenvalue is in fact given by the analytic result. 

Now the limit on the number of bound modes that could be identified in the numerical 
procedure is not really surprising. In the analytic expression the difference between the 
eigenvalues of adjacent modes decreases for the higher modes, whereas in solving the 
eigenvalue problem numerically the differences between eigenvalues that can be resolved 
is dependent on the number of lattice points and the box size. Our use of the non- 
linear tanh map was such that the number of lattice points and box size were increased 
simultaneously. We can see that the eigenvalues of the n = 1, ... ,5 bound modes do 
not change significantly as the number of lattice points is increased. The eigenvalue of 
the n = 6 bound mode, the highest bound mode identifiable with a 100 point lattice, 
does change significantly as the number of lattice points is increased; it tends towards 
the analytic result. The higher the eigenvalue of the bound mode, the further from the 
monopole core the fields of the bound mode extend, hence as the box size is increased 
(note that box size and lattice size increase simultaneously) the number of identifiable 
bound modes increases. We would therefore expect that if we continued to increase the 
number of lattice points (and box size) we would continue to find more bound modes, 
and their eigenvalues would agree with the analytic result. It is worth noting the 1/r 
dependence of the monopole profile H for large r, which in a sense means that all box 
sizes are too small. 

As for the interpretation of the modes, for ipQQ we can interpret the two eigenspectrums 
as being associated with the massless gauge boson and one of the massive gauge bosons. 
The massive gauge boson possesses bound states. For ipQ^ we can interpret the two 
eigenspectrums as being associated with the massless Higgs boson and the other massive 
gauge boson. It is again the massive gauge boson which has bound states. 

To check that these modes are in fact physical we solve the equation for the ghosts. 
For k = the only allowed mode is 

Voo = vT 
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which satisfies the equation 



2(1 - 



This has no bound mode solutions and so the bound states found above are physical. 
This ghost mode cancels the Wi*^ and 05° massless modes. 

These results agree with ref. P where he only studied the ifj^^ perturbations for = 
and not the ip^Q type. He found an infinite number of bound states and numerically 
estimated the lowest to be u"^ ~ 0.768. This agrees well with the value we obtained 
numerically for the first non-zero eigenvalue, but he did not find the bound zero mode. 
This spherically symmetric zero mode is not canceled by the ghosts, and is in fact as- 
sociated with a physical degree of freedom of the monopole. Finally no negative modes 
were found in agreement with the perturbative stability analysis of ref. 



4 Perturbations about the dyon 

Substituting the dyon solution into the gauge fixed second variation of the action gives 



5S, 



9f 



r'^dr I dVt 



{D^)ab + S^' ( ^{H^ - 1) + hA + r-f\(3 - l)H 



and the Fadeev-Popov ghost term as 



(9) 



5S 



FP 



r'dr / cifir/" {D')ab + H^S"' - r^r") 



where the key terms are given by 



ah 



_2 

-V H \ 



w_ 

2" 



- 2^eiabLi + f'f' 



r 



fyc j ^ - - 2) 1 - 6""— 



'W W 



kc 



w 



(A 



H 



r'r IH' {1-W))+ 5*^— (1 - W) 



.H . 



J 



J' - -(1 - w^) + 5''-{i - w) 



0. 



Note that in this case there is a linear time derivative so it is not obvious that it is 
possible to set up a standard eigenvalue problem. 
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The dyon solution is also invariant under rotations generated by and so we again 
expand the perturbations in eigenstates of 



km j=—l 
km j=—l 



E E ^^n^m 
km j=—\ 

2 h 
km ji=0 32=— jl 



and again find that the perturbations separate into two channels under parity 



Ujkm 
km 





\ 


















^rkni 




0fcm 








\ 


J 



\ 0r y 

Substituting these expansions into dSgf and performing the angular integrals gives the 
perturbation operator 



m: 



(±) 



■a/3 



</3 



(±)(2) 
i/3 



+ 2— C, 
r 

H 



+ 4ii(i_iy)cg)(5)+4(i7'--(l-W^))ci^ 



'(±)(4) 
9 

i7 



^(±)(6) 



+ 4^(1 - W)cS»'>e„, + 4 (/ - ^(1 - W)) cS>'«e„, 
for the Higgs and gauge fields, and 



+ 



2W 



(±)(2) 
=/3 
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for the ghosts. We have used the same notation as for the monopole and the non-zero 
elements of the matrices C^^^^~^^ are the same as for the monopole. The non-zero 

elements of the matrices C^^^'^'^^ are the same as those of the matrices Cj^''^^'^^ but with 
the a,P refering to the W fields instead of the fields. We have not calculated the 



elements of for general k, but for the case of interest {k = 0) the elements are all 

zero, and so in this case we can set up a standard eigenvalue problem as for the monopole. 

We will restrict our attention to the Bogomol'nyi limit and consider the analytic 
form of the dyon solution given earlier. Substituting this solution into the perturbation 
equation and rescaling 

C2)-i/\ = (1 + t^{l- C^^'h = (1 + e^'H 

gives Mg^ = (1 + e^y^Jd^^^ where 



2 



(±)(3) ^yV_ ^^^^ _ o^r^(±)(4) 



+ aVTT^ (^^(1 - W)Cl^/'^ + i^H' - ^(1 - W)j cl^p^^'^j 
+ 4e (^^(l - W)di^^\^, + (f - ^{1 - W)) C;$^^%J . 



Note that the quantity we have used in the rescaling (1 — C^) is equal to vT+~e^, 
which is the ratio of dyon mass to the monopole mass. 

k = modes 

If we now consider the k = modes we find that the allowed modes are 



V'oo — ( p^OO ) V'l 



00 



00 



00 



ipQQ satisfies the eigenvalue equation 



where the elements of M are the same as for the monopole, and so this equation is the 
same as equation ^. So we have the same bound state solutions with eigenvalues 

ul = 1 where n = 1, 2, 3, . . . 
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as for the monopole. Note that there is exphcitly no e-dependence of the eigenvalues. 

The equation for ip^Q is different for the dyon and does possess some e-dependence. 
ipQQ satisfies the eigenvalue equation 



where 



00 



.2 , 2(1 - Wf 



M22 = M33 = -V^ + 
M12 = -M21 = e2\/2 ^^^~^^ 



Mi3 - M31 = v^r+^2v^ ^^^ 

r 

M23 = M32 = 
However in the far field limit the perturbation equation becomes 



/ 


_V2 + 1-2 







\ 


1 <° ^ 









-V2 

V ^ 










V 








V J. 


1 


V ) 





which has no e-dependence. So as for the monopole we find that the W^^ field has bound 
solutions with eigenvalues 

oJi — 1 where n = 1, 2, 3, . . . 

The A; = ghost equation for the dyon is the same as for the monopole, and so in the far 
field hmit the ghosts cancel the Wf^ and 0?° massless modes just as for the monopole. 

If we were to solve the eigenvalue equation with the full analytic solution, then 
from our results for the monopole we would expect the eigenvalues to agree with the 
analytic far field result i.e. with no e-dependence. On physical grounds, however, we 
would expect some of the eigenvalue spectra to possess an e-dependence because there 
should be an interaction between the electric charge of the dyon and the electric charge 
of the 1^-boson. The situation is unfortunately obscured by the fact that eigenstates 
of K are not eigenstates of charge. The perturbation equations for the higher k modes 
do in general possess an e-dependence in the far field limit, but we have not solved the 
perturbation equations in these cases. 
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5 Perturbations about the point electric solution 

Substituting the point electric solution into the perturbation equations [l| gives, in the 
BogomoFnyi limit 



and 



where 



[do \ 





±2iDi4 

+ (03)2 ^2iG^ 



Oi 

2 I rJs3\2 



0, 




(10) 



111 



d2 



r 



^-2 



The equations for the gosts {jf ^ 77 ) are 



(d^ + (0^)^) 




0. 



The point electric solution is invariant under the total angular momentum operator J = 
L + S", and so we expand the fields in terms of spherical harmonics Yjm and vector 
spherical harmonics Y^^^ as 



Wo 



^ 4'jmYjm, 
jm 

jm 



jm 



jm s=—l 



where we have suppressed the ±, 3 isospin labels. The angular momentum / is given by 
l=j+s. 
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Unlike the cases of the monopole and dyon, in this case because the total angular 
momentum operator is independent of the isospin generator T we can construct simul- 
taneous eigenstates of J and the charge operator Q = T3. The (0^, Wj^) perturbations 
therefore have zero charge and so the physical modes of equation |10| involve the neutral 
massless gauge boson and the neutral Higgs scalar. The perturbations have 



charges ±1 and so the physical modes of equation |ri| involve the massive charged gauge 
bosons. 



j = modes 

If we consider the j = modes then we find that the only allowed mode is 




which satisfies the perturbation equation 



V 



-(--f)'-v? + (i-f) 

-(^-^Y + 



2iZ/r^ 
- 2iZ/r^ 



-2zZ/r^ -2zZ/r2 - (a; - f ) -V^ + ^ + (l-f)^ 




(12) 



where we have taken the matrix elements for (0+, W^) in equation |Tl|. This is not of the 
form of a standard eigenvalue problem as it stands, but if we rescale p = 2(1 — cu'^Y^'^r 
and define 

(I -co) 



A 



then we obtain 





1 A 
4 p 


2iZ/p^ 


\ 







- + i - 

P ' 4 


A - 2iZ/p^ 




W^oo = 


^ -2iZ/p2 


- 2iZ/p^ 


^ P + 4 p ^ p2 


/ 





(13) 



which is now an eigenvalue problem with eigenvalues A. If we rewrite this as 






( r \ 





















where is the 3x3 identity and the matrix A is given by 



A 









-z 








z 


-z 


-z 


-1 
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then the problem of diagonahzing the perturbation matrix has been reduced to diago- 
nahzing the matrix which is easy. The eigensolutions of 

A e = Q e 



are 



Qi = ei 



/ 1 
1 



= 62 



/ 1 
-1 

V 



= -1 63 



/ z 

-z 

\ 1 



The corresponding eigenvectors of the perturbation matrix are 









( ^ ) 




( ^ \ 


#00 






-1 




-z 








I ) 




[ 1 / 



where the functions Wi{p) satisfy 



for i = 1,2 



If we write Wi{p) — Ui{p)/p, then the equations for the Ui are of the form of Whittakers 
equation 

fd^ 1 



dp^ 4 ^ ^ 



A l-^i'^ 



Ui{p) = 0. 



p p- 

where p — 1/2 for ui and 1*2, and p — 3/2 for u^. As we stated earher, for a;^ > 1 
(imaginary p) these equations have Whittaker function solutions Wts^^^j,{p) for all A which 
are oscillatory at large p. For 1 > o;^ > (real p) the Whittaker fuction solutions only 
exist for 



h. — rir + p-\- 



0,1,2,3,... 



and the solutions decay exponentially at large p. 

At this point it is worth considering the ghosts. The corresponding equation for the 
J = mode of the charged ghosts is 



^ Z^ 2' 
r 



%o = 0. 



(14) 



If we similarly rescale p = 2(1 — uPY^'^r and define 



A 
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then we can rewrite the ghost equation as 



which is the same as the equation for the wi and W2 eigenmodes. The ghosts therefore 
cancel the wi and W2 eigenmodes leaving the W3 eigenmode as the physical mode. 
For the W3 mode we have = 3/2 and so there are physical bound modes for 

A = n = 2,3,4, ... 

Substituting this into the expression 



A = Z 



2 



and squaring it gives a quadratic equation in uj with the two solutions 

(1 — In^) 

a;+ = 1, and = - n = 2,3,4, ... 

(1 + jn^) 

For our point electric solution to correspond to a boson we take Z = l/47r in our 
rescaled units. Since then = l/lGir"^ < 1 we can approximate a;_ as 

/ 2zn 

^---(1-^1 n = 2,3,4,... 
Now the perturbations are eigenstates of charge and so we have 

So the solution therefore corresponds to the boson and the U-. solution corre- 
sponds to the W~ boson. Our result therefore agrees with the result of ref. p2| where it 
is shown that the force between two gauge bosons of the same charge is zero, and hence 
no bound states, whilst the force between gauge bosons of opposite charges is twice the 
electrostatic coulomb force. 

We note that the point electric solution used here has no perturbative instabilities 
[uj imaginary), no matter how large the charge AttZ. This differs from the point electric 
solution 

= --, <j)^ = r] 
r 



which possesses perturbative instabilities above a critical charge |23 
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6 Discussion 



We have found the eigenvalue spectra for spherically symmetric perturbations about the 
monopole, dyon and point electric solutions. For the perturbative modes of positive parity 
about both the monopole and dyon we argue that the eigenvalue spectra of bound modes 
are the same using both the far field form and the full form of the Prasad- Sommerfield 
solution. These bound states are interpreted as bound states of the massive H^-bosons 
in the presence of the monopole and the dyon. The force between the monopole and 
the massive gauge bosons consists of two components. One part is the electromagnetic 
interaction between the monopole magnetic charge and the magnetic dipole of the gauge 
boson. However it was shown in ref. |^ that this interaction by itself is not strong enough 
to give bound states. The other part of the monopole-gauge boson force arises from the 
Higgs field of the monopole. Because the gauge boson acquires its mass from the Higgs 
field expectation value, the gauge boson has a lower energy when it is in the presence 
of the monopole core and this gives rise to a force which is responsible for confining 
the gauge boson to the core of the monopole. The force between the dyon and the 
massive gauge bosons has an electrostatic force in addition to the monopole-gauge boson 
force. The effect of this electrostatic force is unfortunately obscured by the fact that 
the perturbative modes about the dyon (and monopole) are not in eigenstates of charge. 
This is not due to any failings of our perturbative method but is an intrinsic property of 
both the dyon and monopole solutions. 

The perturbations about the point electric solution are in eigenstates of charge. The 
force between two charged gauge bosons consists of an attractive Higgs force and an 
electrostatic force of equal magnitude. For gauge bosons of equal charge the two forces 
cancel and in this case we find no bound states, whereas for oppositely charged bosons 
the two forces add and our results for the bound states in this case are consistent with 
this. 

We set out to answer the question: are the bound state energies of the various two 
particle systems of the Yang-Mills-Higgs model related by electromagnetic duality? We 
are specifically considering the SL{2,Z) symmetry of Sen [jl3[. The original conjecture 



of Montonen and Olive |TD[ was that there are two equivalent choices of action, either 
with the ly-bosons as the gauge particles and the monopoles as the sohtons, or with the 
monopoles as the gauge particles and the VT-bosons as the solitons. The spectrum of 
dyons does not, however, fit nicely into this picture. In the electromagnetic duality of 
Sen Jl^ the electric (q) and magnetic (g) charges of a particle state are given by 



q + ig = qoirriT + n) m,n & Z 
where go is the minimum electric charge and 

e i 

6' is a new parameter of the theory. Plotting q + ig in the complex plane gives the charge 
lattice {q,g), and the single particle states correspond to the primitive vectors of this 
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charge lattice. A basis for the charge lattice is given by any two non-collinear primitive 
vectors, and choices of basis correspond to choices of equivalent actions. Invariance of 
the theory under electromagnetic duality is now realised as invariance of the theory to 
changes of basis of the charge lattice. 

If we choose the primitive vector go to correspond to the electric gauge boson and the 
primitive vector q^r to correspond to the magnetic monopole, then a change of basis to 
(?o, I'o^') is given by fH 






where a,b,c,d ^ Z and ad — be = 1, so the matrices A form the group SL{2, Z). Under 
this transformation r becomes 

, ar + b 

^ = 1- 

CT + a 

We now consider whether our results are consistent with this picture. Our H^-boson- 
monopole system corresponds to choosing the primitive vectors qo as the IV-boson and 
QoT = i/lo SiS the monopole. In this case we find that the spherically symmetric bound 
perturbative modes about the monopole have energies 



11? J \ m? 



where q is the charge of the ly-boson and g is the charge of the monopole. If we now 
apply the SL{2, Z) transformation 






then the ly-boson still corresponds to the point qo, but the monopole in the new basis 
corresponds to the point q^r + q^, which is the dyon of the old basis. The product of 
the M^-boson charge q and the monopole charge g is invariant under this transformation, 
so the bound state energies of the two particle system should be invariant under this 
transformation. We did indeed find that the VT-boson-dyon bound state energies are 
given by the same expression, despite our reservations about the VT-field perturbations 
not being in eigenstates of charge. 
The SL{2, Z) transformation 




r 
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applied to the original ly-boson-monopole basis {qo,i/qo) corresponds to the duality 
transformation of Montonen and Olive [|10|- The transformation displays the interchange 
of weak and strong coupling r —'^/t- This weak-strong coupling interchange also 
occurs for the transformation 




t' 





where the new gauge boson again corresponds to the old monopole but the new monopole 
corresponds to the old dyon with charges (—go, i/lo) in the original basis. The product qg 
is again invariant under this transformation and so we would expect the monopole-dyon 
bound state energies to again be given by 



The monopole-dyon two particle system was studied in ref. |]14| where they found physical 
bound states for 

1 

' = 2' 



^^2_^2^1/2. 2_^2^1/2^ 



4^2 



for > 



where s, n are integers and s is the relative electric charge between the monopole and 
dyon. This is obviously different from what we have above. 

Both our results and those of ref. are only approximations, and under the SL{2, Z) 
duality transformations they correspond to different limits. Our approach was perturba- 
tive, and we found the harmonic oscillator states about the monopole and dyon. Although 
in the quantum theory these quantized harmonic oscillator states are interpreted as the 
gauge bosons of the theory, this approach only accurately describes the VT-boson-soliton 
two particle system in the limit <^ 1. 

In ref. [^] they considered the dynamics of two non-relativistic BPS monopoles us- 



ing the Atiyah-Hitchin metric on the space of collective co-ordinates of the monopoles. 
These collective co-ordinates consist of the three spatial positions and a phase. The 
spherically symmetric zero mode we obtained for both the monopole and dyon corre- 
sponds to changes in this phase. In the limit of large monopole separation they used 
the Taub-NUT metric as an approximation to the Atiyah-Hitchin metric, and solved the 
non-relativistic Schrodinger equation in Taub-NUT space for the relative motion of the 
two monopoles. 



The result of ref. \14] is non-relativistic and was obtained in the limit of weak electric 



coupling. If we consider the energy of the dyon given earlier in the limit g -C 1 we find 
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So the total energy of the static monopole-dyon system in the hmit of infinite separation 
for -C 1 is 

----- 

If we exphcitly include factors of h so that the gauge coupling is hq, then the weak 
coupling condition is hq <^ 1. Now we have taken ^ = 1, and so weak coupling condition 
is g <^ 1, whereas in ref. ||14| they use units q = 1 and so the weak coupling condition is 



given hy h <^ 1. With the units g = 1, = 1 of ref. the energy given above becomes 

E - 




and in the centre of motion frame with the mass and length in units of An, bound states 
occur for 

< e < 

4 

Now the electric charge of a dyon in the presence of a monopole is not a good quantum 
number, but the relative electric charge (s) is. If, in the centre of motion frame, we 
assign electric charges {Qi = —s/2) to the monopole and {Q2 = s/2) to the dyon, then 
the result 



- 1- 



is consistent with an electrostatic interaction of the form Q1Q2/T between the two par- 



ticles p5 |. 



So both here and in ref. |jT^ bound states due to a Coulomb interaction are found in 
the limit q^ <^ 1. However, the SL{2, Z) transformation between the M^-boson-monopole 
and monopole-dyon systems involves a weak-strong coupling interchange. So our result 
for the M^-boson-monopole system, which holds for ^ 1, would correspond to the limit 
g^ 3> 1 for the monopole-dyon system, and the result of ref. [Q, which holds for g^ <C 1 
for the non-relativistic monopole-dyon system, would correspond to the limit g^ ^ 1 for 
the non-relativistic ly-boson-monopole system. 



In ref. [25| Manton showed that the monopole-dyon interaction could be understood 
in terms of the motion of classical point particles with electric and magnetic charges. If 
we consider a dyon with electric charge g and magnetic charge then the scalar charge 
must be {g'^ + g^)^/^ [^]. So the Coulomb interaction between two static particles with 



charges (gi,5'i) and (g2,5'2) is (HI 

V 



9192 + qiq2 - {9I + qlY'\9l + qIY'^ -■ 

J f 

Now if we consider a monopole (0, g) and a IV-boson (g, 0) we find 
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and for a dyon (g, g) and a VT-boson (g, 0) we find 



~ roi 

r 

which is the same as the VT-boson-monopole interaction for <^ 1 (note g = 1/q). The 
monopole-dyon interaction is given by 

1 

r 



V 



,2n1/2 



9 -gig +q 

2 



q^ 



for 

2r 



which gives the result for n ^ s of ref. |jT^. For q ^ g we find 



V - 



ig 



which is the same as for the VT-boson-monopole system. This simple argument supports 
the view that our result for the bound state energy 



1 - 



Qg 



applies for g <^ 1 for the relativistic 14^-boson-monopole and W^-boson-dyon systems, and 
also for the relativistic monopole-dyon system for g^ ^ 1. 

If the bound state energies of these three two-particle systems are indeed related 
by SL{2, Z) transformations on the charge lattice, then the result is in fact far more 
general. This is because any two non-coUinear primitive vectors of the charge lattice 
can be transformed to the basis {qo,i/qo)- Our result then gives the relativistic bound 
state energies in the weak coupling limit (g^ <^ 1), and the result of ref. jMl gives the 
non-relativistic bound state energies in the strong coupling limit (g^ ^ 1). 

If we consider two general single particle states with the same charges {g,q), they 
exert no static Coulomb forces on each other since 



V 



2 I 2 

g +q 



{g' + q'f'\g' + q 



2^l/2 



1 

- = o, 

r 



and so they form no bound states with one another. Two general single particle states 
with opposite charges, however, have a static Coulomb force 



V 



-g -q 



{g' + q'f'\g' + q'f' 



-2 



g^ + q^ 



and can form bound states. By an SL{2, Z) transformation we can transform this two 
particle system to the W~^-W^ boson system, which we found to have bound state 
energies 

(1 - gVn^) 



UJ 



(l + gV«') 



for g ^ 1. 
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As to our original question of whether the bound states of the various two particle 
systems are related by electromagnetic duality, we have no definitive answer but our 
results would appear to be consistent with it. 

Finally we note the possible existence of an unexpected coincidence. Our numerical 
solution of the perturbations about the monopole gave the lowest massive bound mode to 
have an eigenvalue of iJ^ ~ 0.768 vn? in the Bogomol'nyi limit. We, however, argued that 
the omission of the origin in our numerical method is responsible for it being different 
from the analytical far field result of cu^ = 0.75 vr? . In ref. [0 we numerically solved for 
the perturbations about the local string, and in the Bogomol'nyi limit imu = ttlw) we 
found a bound boson with eigenvalue a;^ ~ 0.775 m^. Our numerical method in this case 
also necessitates the omission of th origin due to a centrifugal term. In addition 
the string profiles had to be solved for numerically because there is no analytic solution. 
The fields of this bound mode vary in exactly the same place that the string profiles 
vary most, and so the value of this eigenvalue is exceptionally sensitive to the accuracy 
of the string profiles. So we argue that our numerical result of uj^ ~ 0.775 would be 
consistent with the true eigenvalue being u'^ = 0.75 m^. If this is the case then the kink, 
string and monopole all possess bound bosons with an eigenvalue of cj^ = 0.75 when 
they are in their respective Bogomol'nyi limits [ill {f3 = rn^j/m^ for the monopole; 
13 = m|^/m^ —>■ 1 for the local string; and all values of the Higgs self coupling for the 
kink). It is precisely when these theories are in their respective Bogomol'nyi limits that 
supersymmetric extensions can be constructed where the Bogomol'nyi bounds |Tl|] can be 
derived algebraically from the supersymmetry algebra. So it may be that this apparently 
common bound of tu^ > 0.75 for the bosons of the theory in the presence of a defect 
with topological charge one, could be similarly derived in the supersymmetric extensions 
of the respective theories. 
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Table captions 

Table 1: Diagonal elements of Cl^p^^~^^ 

Table 2: Off-diagonal elements of C^g'"'"'^'* 

Table 3: Off-diagonal elements of C^^'"-^"^'* 

Table 4: Elements of C^f ^^"^^ 

Table 5: Bound mode eigenvalues for the monopole 
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